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PERIODIC-COEFFICIENT DAMPING ESTIMATES, AND STABILITY OF 
LARGE-AMPLITUDE ROLL WAVES IN INCLINED THIN FILM FLOW 

L.MIGUEL RODRIGUES AND KEVIN ZUMBRUN 


Abstract. A technical obstruction preventing the conclusion of nonlinear stability of large-Froude 
number roll waves of the St. Venant equations for inclined thin film flow is the ’’slope condition” of 
Johnson-Noble-Zumbrun, used to obtain pointwise symmetrizability of the linearized equations and 
thereby high-frequency resolvent bounds and a crucial nonlinear damping estimate. Numerically, 
this condition is seen to hold for Eroude numbers 2 < F ^ 3.5, but to fail for 3.5 ^ F. As hydraulic 
engineering applications typically involve Froude number 3 ^ E ^ 5, this issue is indeed relevant to 
practical considerations. Here, we show that the pointwise slope condition can be replaced by an 
averaged version which holds always, thereby completing the nonlinear theory in the large-E case. 
The analysis has potentially larger interest as an extension to the periodic case of a type of weighted 
“Kawashima-type” damping estimate introduced in the asymptotically-constant coefficient case for 
the study of stability of large-amplitude viscous shock waves. 


1. Introduction 

The St. Venant equations of inclined thin film flow, in nondimensional Lagrangian form, are 

dtT - dxU = 0 , 

^ ^ + =l-Tu^+ ud^{T~‘^dxu), 

where r = l//i is the reciprocal of fluid height h, u is tangential fluid velocity averaged with 
respect to height, x is a Lagrangian marker, E is a Froude number given by the ratio between a 
chosen reference speed of the fluid and speed of gravity waves, and v = with Re the Reynolds 
number of the fluid. The terms 1 and tu^ on the righthand side of the second equation model, 
respectively, gravitational force and turbulent friction along the bottom. Roughly speaking, F 
measures inclination, with F = 0 corresponding to horizontal and E ^ oo to vertical inclination of 
the plane. 

An interesting and much-studied phenomenon in thin film flow is the appearance of roll-waves, 
or spatially periodic traveling-waves corresponding to solutions 

(1.2) {T,u){x,t) = {f,u){x — ct) 

of (1.1). These are well-known hydrodynamic instabilities, arising for (1.1) in the region E > 2 
for which constant solutions, corresponding to parallel flow, are unstable, with applications to 
landslides, river and spillway flow, and topography of sand dunes and sea beds [BM04]. 

Nonlinear stability of roll-waves themselves has been a long-standing open problem. However, 
this problem has recently been mostly solved in a series of works by the authors together with 
Barker, Johnson, and Noble; see [JZNll, BJRZll, BJN+10, JNRZ14, BJN'’'15]. More precisely, it 
has been shown that, under a certain technical condition having to do with the slope of the traveling- 
wave profile (t,u), spectral stability in the sense of Schneider [Sch98, Sch96, JZNll, JNRZ14], 
implies linear and nonlinear modulational stability with optimal rates of decay, and, moreover. 
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asymptotic behavior is well-described by a system of second-order Whitham equations obtained by 
formal WKB expansion. 

In turn, spectral stability has been characterized analytically in the weakly unstable limit F —)■ 2 
and numerically for intermediate to large F in terms of two simple power-law descriptions, in the 
small- and large-F regimes, respectively, of the band of periods X for which roll waves are spectrally 
stable, as functions of F and discharge rate q (an invariant of the flow describing the flux of fluid 
through a given reference point) [BJN+lh]. That is, apart from the technical slope condition, there 
is at this point a rather complete theory of spectral, linear, and nonlinear stability of roll wave 
solutions of the St. Venant equations. However, up to now it was not clear whether failure of the 
slope condition was a purely technical issue or might be an additional mechanism for instability. 

Precisely, this slope condition reads, in the Lagrangian formulation (1.1)-(1.2), as 

(1.3) 2uu^ < F-\ 

where u is the velocity component of traveling wave (1.2). It is seen numerically to be satisfied 
for F ^ 3.5, but to fail for F ^ 3.5 [BJN'’'15]. For comparison, hydraulic engineering applications 
typically involve Froude numbers 2.5 ^ F ^ 20 [AeM91, Bro69, Bro70]; hence (1.3) is a real physical 
restriction. From the mathematical point of view, the distinction is between small-amplitude, slowly 
varying waves for which (1.3) is evidently satisfied and large-amplitude, rapidly-varying waves, such 
as appear in the spectrally stable regime for small and large F, respectively [BJN+lb]. 

The role of condition (1.3) in the stability analysis is to obtain pointwise symmetrizability of the 
linearized equations and thereby high-frequency resolvent bounds and a crucial nonlinear damping 
estimate used to control higher derivatives in a nonlinear iteration scheme. The purpose of the 
present brief note is to show, by a refined version of the energy estimates of [JZNll, BJRZll], 
that the pointwise condition (1.3) can be replaced by an averaged version that is always satisfied, 
while still retaining the high-frequency resolvent and nonlinear damping estimates needed for the 
nonlinear analysis of [JZNll, JNRZ14], thus effectively completing the nonlinear stability theory. 

The remainder of this paper is devoted to establishing the requisite weighted energy estimates, 
first, in Sections 2-3, in the simplest, linear time-evolution setting then, in Sections 4.1 and 4.2, 
respectively, in the closely related high-frequency resolvent and nonlinear time-evolution settings. 
The estimates so derived may be seen to be periodic-coefficient analogs of weighted “Kawashima- 
type” estimates derived in the asymptotically-constant coefficient case for the study of stability 
of large-amplitude viscous shock waves [Zum04, Zum07, GMWZ06], to our knowledge the first 
examples of such estimates specialized to the periodic setting. We discuss this connection in 
Sections 5 &: 6. More, this seems to be the first instance of a proof of hypocoercive^ decay where 
periodicity is used in a crucial way. We note, finally, the relation between these weights and the 
“gauge functions” used for similar purposes in short-time (i.e., well-posedness) dispersive theory 
[LP02, BGDD06, Miel5], a connection brought out further by our choice of notation in the proof. 
This indicates perhaps a potential for wider applications of these ideas in the study of periodic 
wave trains. 


( 2 . 1 ) 


2. Preliminary observations 

Making the change of variables x —>■ x — ct to co-moving coordinates, we convert (1.1) to 

Tt - CTx -Ux = 0, 

Ut - CUx -b ((2F^)“V“^)2, = l-TU^ + 


^The reader interested in replacing Kawashima-type estimates in the more general context of hypocoercive decay 
estimates is referred to [Vil09, Remark 17] and references in [Rodl3, Appendix A], especially [BZll]. 
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and the traveling-wave solution to a stationary solution U{x,t) = (t(x, t),tt(x,t)) = ('f(x),n(x)) 
convenient for stability analyis. 

We note for later that the traveling-wave ODE becomes 

( 2 . 2 ) -cfx-Ux = 0, -cUx + {{2F‘^)~^f~‘^)x = l-fu‘^ + i^{f~‘^Ux)x, 

yielding the key fact that 

(2.3) f{f)ux = cf{f)fx 

is a perfect derivative for any function /(•), hence zero mean over one period. We note also as in 
[JZNll] that c 7 ^ 0, else u = constant and the equation for r reduces to first order, hence does not 
admit nontrivial periodic solutions. Linearizing about U = (r, u) gives the linearized equations 

Tt - CTx -Ux = 0, 

(^•^) / \ (—2 \ -2 9-- 

Ut — CUx — {CtTjx = V[T Ux)x — U T — lUTU, 

where 


(2.5) a := T ^{F ‘^ + 2i'Ux)- 

With this notation, the slope condition of [JZNll] appears as f^a > 0. We note that, by (2.3), 
the mean over one period of g{f)a is positive for any positive g: 

(2.6) {g{f)a) = {g{f)f~^F~^) > 0. 


That is, (any reasonable version of) the slope condition holds always in an averaged sense.^ An 
approximate asymptotic diagonalization in the large spectrum regime- see [BJRZll, BJN+IS], in 
particular [BJN+15, Appendix A]- reveals that the sharp^ relevant averaged conditions is 


(ar^) 


> 0. 


V 

We shall show in the rest of the paper that this averaged condition is in fact sufficient for the 
nonlinear analysis of [JZNll, JNRZ14]. 


3. Linear damping estimate 


Introduce now some ‘gauge’ functions (/>i, (j )2 and (ps and define for U = {t,u) the energy 
(3.1) £{U) := j + + (psTUx^ . 

A brief computation yields that solutions U of (2.4) satisfy 

2 ) ^ ~ + "'^ 3 ) rl - {^(p2) ul^ + {(pi - a<p2 + ^(ps) TxUxx) 

+ 0 {{\\ u \\ jj2 + ||r||jyi)(||tt||j^i J- ||r||i2)). 


The original gaugeless strategy that works when a is positive may be achieved by choosing (pi = 1, 
(p2 = and 0 < (/>3 = constant <C 1. The possibility of choosing (p2 = (pi /a while keeping 

both (pi and (p 2 positive is a direct manifestation of the fact that in this case the first-order part of 
system (2.4) is symmetrizable. For the general case, of interest here, we instead take 

(3.3) £(*), . 


-)')>i - 0. .^ 1 ( 0 ) > 0, 


^Here and elsewhere we use {h) to denote mean over one period of a function h. 

^In the sense that there exist curves of spectrum for the operator L- defined below- that are going to infinity and 
whose real parts converge to —{aF)lv. 
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(3.4) 


(pi — q ;(/>2 + ^(/>3 = 0, 0 < (/>2 = constant <C 1, 

so that (/>3 is chosen to kill the indefinite cross-term and the fact that (pi is not constant and thus 

does not commute with the generator of system (2.4) is used to average and cancel the “bad” 
_2 

oscillating part of through the arising nontrivial commutator. 

With these choices, we obtain after another brief computation 


(3.5) 





</>2 rj + 


t4>2) 


+ 0{\\u\\h^ -I- ||'r||j^i)(||u||j^i -I- ||r||£,2)) 

< —11i{\\Uxx\\‘l2 + || Ta ;||| 2 ) -I- Cl (||u||j^2 -|- ||r||j^l 


1^11^1 



II1 /2 11 111 /2 

for some positive rji and Ci, whence, by interpolation inequality ||u||hi < \\u\\j ^2 IfII /^2 and the 
fact that S{U) ~ (||ua ;||^2 + ||' 5 e||^ 2 ) modulo ||t||^ 2 , 

( 3 . 6 ) ^ -v£mt)) + cwumh, 

for some positive tj and C, a standard linear damping estimate. 

Note that, in the step £{U) ~ (||ua ;||£2 + ||tx||^ 2 ) modulo L^, we have used in a critical way that 

C/af^ {af‘^)\ 

J \ u V )’ 

hence (pi and l/<^i, remains bounded, a consequence of periodicity plus zero mean. 

To derive a first corollary from the key estimate (3.6), one may combine it with the standard 
bound 

( 5 ^^ + 5 “^) = - J + 0{{\\u\\hi+ \\t\\l2){\\u\\l2+ \\t\\l2)) 

to obtain the following lemma. 


Lemma 3.1. There exist positive 6 and C such that any U solving (2.4) satisfies for any t > 0 

||C/(t)||Hi < Ce-^^\\Um\Hi + C fe-^^^-^^\\U{s)h2ds. 

Jo 


4. Applications 

Lemma 3.1 by itself is not of much direct practical use. However, as we will now show, we 
can readily adapt its proof, and especially estimate (3.6), to obtain various useful forms of high- 
frequency damping estimates. The reader unfamiliar with these considerations may benefit from 
first having a look at [Rodl3, Appendix A] for a terse introduction to this approach. Indeed what 
follows stems directly from the mere introduction in the classical strategy described there of gauges 
leading to (3.6). See also [Zum04, Zum07] for related estimates in the shock wave case. 

4.1. High-frequency resolvent bound. An important part of the proofs in [JZNll, JNRZ14] is 
dedicated to estimates of semigroups generated by linearization around a given wave, to be used 
in an integral formulation of the original nonlinear systems. These estimates are deduced from 
spectral considerations and the noncritical part of the linearized evolution is directly controlled by 
an abstract spectral gap argument that only requires uniform bounds on certain resolvents. Our 
claim is that a spectral version of Lemma 3.1 does provide these uniform bounds. 

To be more specific let L denote the operator generating the linearized evolution around U, that 
is, such that system (2.4) reads Ut — LU. The operator L is a differential operator with periodic 
coefficients but acting on functions defined on the full line. We do not apply directly spectral 
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considerations to L but rather to its operator-valued Bloch symbols L^, associated with the Floquet- 
Bloch transform- see [JNRZ14, Rodl3] for instance. Explicitly, if H denotes the fundamental 
period of tj, for any Floquet exponent ^ in the Brillouin zone [—tt/S, vr/S), the operator acts 
on functions of period S by := 

The operator dx itself has Bloch symbols dx + As a result, when dealing with L^, the 

(equivalent) norm of interest on F^g^(0,S) is ||/||h| = (Z]fc=o IK^a; + ^0^/lli2(o =)) ^ • Consider 
now the resolvent equation 

(4.1) {\-L^)U = F. 

Letting (•, •) denote complex inner product, we find by computations essentially identical to those 
in Section 3, substituting XU for Ut and dx + for dx, that, defining (f>j as in (3.3)-(3.4), and 

£^{U) := ^{<j)i{dx + iC)T, {dx + i^T) + ^{(^)2{dx + i^x + iC)u) + {dx + iC)u), 

one derives 

2iR(A) £^{U) = ^{2X£^{U)) 

= ^{{(j)i{dx + iO'u, {dx + iO{^'>J-)) + {4>2{dx + iCh, {dx + iC){^T)) 

+ {4>3T, {dx + + {4>3{dx + iOu, Ar)) 

^^ 2 ) = - j{^{4>i)x + a(f>s)\{dx + i0^f - J{^(f>2)\{dx + i^fu\‘^ 

+ "*^2 + [{dx + ii)T {dx + 

+ 0(||m||j^ 2 -I- ||r||j^i)(||u||j^i -I- ||r||i2)) + 0(Ill’ll//! IIE||//i) 

<-r?£:^(t/) + C(||C/||i. + ||Fl|^p. 

for some positive r] and C uniform with respect to A and 
Combining (4.2) with the easy estimate 

mh = ^ |(C/,A[/)|<C|A|-i(||[/||^i + ||F||i.), 

for some C, we obtain for > —ri/2 and |A| sufficiently large, the estimate 

\mi, < c\\F\\i, 

for some uniform C. Incidentally, since has compact resolvents hence discrete spectrum com¬ 
posed entirely of eigenvalues, the above estimate also implies that such A do not belong to the 
spectrum of L^. 

More generally, by adapting the previous computations to higher-order estimates, along the lines 
of the method expounded in next subsection, one proves the following result required by the analysis 
of [JZNll, JNRZ14]. 

Proposition 4.1 (Resolvent bounds). For any positive integer s, there exist positive p, C and R 
sueh that if X G C is such that |A| > R and 3?(A) > —rj then, for any [—vr/S, vr/S), we have 

^ ^ ^Hi,^^{o,'B.){Li) and ||(A - L^)“^ < C. 

This offers a direct replacement for [JZNll, Appendix B] without assuming any condition on the 
background wave U. 
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4.2. Nonlinear damping estimate. The other place where high-frequency estimates play a role 
in the arguments of [JZNll, JNRZ14] is in providing a nonlinear slaving bound that shows that 
high-regularity norms are controlled by low-regularity ones and enables us to close in regularity a 
nonlinear iteration. With the strategy implemented above we are also able to reproduce this bound 
without assuming the slope condition (1.3). 

To be more specific let us first warn the reader that, because of the complex spatio-temporal 
dynamics that take place around periodic waves, the appropriate notion of stability is neither the 
standard one nor the simpler orbital stability but space-modulated stability, as recalled in the next 
subsection. For this reason, following [JZNll, JNRZ14], instead of directly estimating U — U, where 
U = {f,u) is a solution of (2.1), we need to introduce (R,V') such that 

(4.3) V{x,t) = U{x —'ijj{x,t),t) — U{x), 


intending to prove that V and the derivatives of remain small provided that they are sufficiently 
small initially. Mark that even if ^ is initially zero, as assumed in [JZNll], one may not achieve 
the latter goal while imposing ip = 0. In other words a modulation in space, encoded by a space- 
time dependent phase is in any case needed in the argument. See the detailed discussions in 
[JNRZ14, Rodl3]. Our new unknowns, which have to be determined together in a nonlinear 
way, are then V = (r, u) and ip, and a specific educated choice, that we shall not detail here, 
is then needed to obtain concrete equations for those. However, let us at least mention that in 
constructions of [JZNll, JNRZ14] the phase shift ip is always slow so that only high regularity 
control on V remains to be proved. This is what we provide now. 

To do so in a precise but concise way, we set /(r) = (2F^)“^r“^, g{T) = vt~‘^ and h{T,u) = 
1 — Tu^. Then U in (4.3) solves (2.1) provided that V = (t,u) and ip satisfy 
(4.4) 


{l-iP^)Vt-LV = 


+ 


1pt{U + V)^ + + (^_^g(^T + T){u + U),^ ^ 

( ° ^ 

\{g{T + t)- g{T))u,^ -I- {g{T + t) - g{T) - g'{T)T)u,j; J ^ 

(/(r + r)-/(f)-/'(r)r)^ + {h{U + U) - h{U) - dhiU){U)) 


Defining the modified energy 


(4.5) := ip^) (^\(PitI + \(p2T‘^ul + epzTUa^ , 

repeating the argument of Section 3, absorbing nonlinear terms into the linear ones and separating 
out ip terms using Sobolev’s embeddings in Gagliardo-Nirenberg form and Young’s inequality, we 
obtain, in analogy to (3.6), that solutions to (4.4) satisfy the nonlinear estimate 


j^£i>{V)<-g£^{V) + C{\\V\\l, + \\iiPt,iPMH^), 

for some positive C and g, provided that we know in advance some sufficiently small upper bound 
on ||(R, ipt,ipx)\\H^ thus are allowed to use Lipschitz bounds for /, g and h and their derivatives 
on a fixed neighborhood of U. Differentiating the equations and performing the same estimate on 
d^V, with higher-order interpolation inequalities, we obtain likewise when k is a positive integer 

(4.6) < -v£Mv) + + ll(V’i,V’a.)||^.) , 

so long as ||R||//i and IKV’x, V’t)!!//'= remain sufficiently small. 

6 



Applying Gronwall’s inequality and recalling that £^{d^V) \\dlV\\l, modulo lower-order 

terms, with constants uniform with respect to ipx satisfying constraints above, we obtain the follow¬ 
ing key estimate showing that higher Sobolev norms are slaved to ||1 ^||l 2 and ||(V't, 

the final nonlinear estimate needed for the analysis of [JZNll, JNRZ14]. This provides a result anal¬ 
ogous to [JNRZ14, Proposition 2.5] and directly replacing [JZNll, Appendix A], without making 
any use of a pointwise symmetrization hence dropping the slope constraint (1.3). 

Proposition 4.2 (Nonlinear damping). For any positive integer s there exist positive constants 6, 
C and £ such that ifV and V' solve (4.4) on [0, T] for some T > 0 and 

sup ^ 

t&[0,T] 

then, for all 0 < t < T, 

(4.7) Ib(t)llHRR) ^ Ce-^^\\v{0)\\jjs(^^) + C j (ll^'(s)|li2(R) + IKV'i, V'x)(s)|li.(R)) ds. 

4.3. Asymptotic stability. As discussed with great detail in [JNRZ14, Appendix D], uniform 
resolvent bounds of Proposition 4.1 and nonlinear slaving estimates of Proposition 4.2 are the only 
structural conditions needed to apply almost word-by-word the arguments of [JNRZ14] to a periodic 
wave of a given ’parabolic’ system. Our foregoing analysis shows that system (1.1) satisfies those 
around any given wave so that all conclusions of [JNRZ14] apply to any spectrally-stable periodic 
wave of (1.1). In particular, any spectrally-stable roll-wave is also nonlinearly-stable, provided that 
one uses definitions of stability adapted to periodic waves of parabolic systems, as we now briefly 
recall. 

A given periodic wave solution to (1.1) N, of period S, is said to be diffusively spectrally stable 
provided that the generator L of the linearized evolution and its Bloch symbols L^, as defined in 
Subsection 4.1, satisfy 

(Dl) a{L) c {A I < 0} U {0}. 

(D2) There exists 6 > 0 such that for all ^ G [—7r/H,7r/H) we have cr{L^) C {A | J?A < — 

(D3) A = 0 is an eigenvalue of Lq with generalized eigenspace of dimension 2. 

(H) With respect to the Floquet exponent derivatives at 0 of the two spectral curves passing 
through zero are distinct. 

From the pioneering work [Sch98, Sch96] to the recent [JNRZ14], conditions (D1)-(D3) have slowly 
emerged as essentially sharp spectral stability conditions for periodic waves of dissipative systems. 
Some form of (H) is also needed but the present form could well be slightly relaxed in a near future, 
see precise discussion in [Rodl3, Chapter 5]. All together, conditions (D1)-(D3) and (H) express 
that the spectrum of L is as noncritical and nondegenerate insofar as allowed by the presence 
around 17 of a two-dimensional family of periodic waves. 

The spatial complexity of the periodic background U precludes any hope for a simple notion 
of nonlinear stability. Over the years there has arisen the concrete remedy implemented in (4.3), 
consisting in introducing a space-time dependent phase shift, though with various possible strategies 
in the prescription of separate- but coupled- equations for V and if. One obvious inspiration for 
introducing a phase in the nonlinear study comes from classical analysis of simpler, asymptotically- 
constant patterns such as fronts, kinks, solitary waves or shock waves, for which the relevant notion 
of stability- orbital stability- already requires the introduction of a time-dependent phase. As 
formalized in [JNRZ14] the corresponding notion of stability for periodic waves- space-modulated 
stability- is obtained by measuring proximity of a function u from a function u in a given functional 
space X with 

6x{u,v) = inf ||uoT-u||x -k ||52,(T - Id)||x. 

'5 
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and not with ||ii — z;||x- At a given time this allows for a space-dependent phase synchronization 
provided that the synchronization differs from the identity by a sufficiently slow phase shift. The 
interested reader is again referred to [JNRZ14, Rodl3] for a detailed discussion of this concept. 
However, we stress again here that there is no hope for a better notion of stability unless the original 
system exhibits some nongeneric null conditions, denoted “phase uncoupling” in [JNRZ14]. 

With these definitions in hands, our analysis combined with the arguments of [JNRZ14] yield 
the following stability result. 

Theorem 4.3 (Nonlinear stability). For any integer K, K > A, a diffusively spectrally stable 
periodic wave of (1.1) is nonlinearly asymptotically stable from L^(M) n i4^(M) to H^{W) in a 
space-modulated sense. 

More explicitly, ifU satisfies (D1)-(D3) and (H) then, for any K > A, there exist positive e and 
C such that any Uq such that 5iif^HK{Uo, U) < e generates a global solution U to (1.1) such that 

Vt€M+, 5HK{U{;t),U)< C 5LinHK{Uo,U) 

and 

6HK{U{;t),U) 0 . 

The actual proof provides a much more precise statement including, for instance, a bound of 

dLp{U{-,t), U) by C (1 -|-similar to those for L^-norms of a heat kernel or of self-similar 
solutions of viscous Burgers’ equations. The reader is referred to [JNRZ14, Theorem 1.10] for 
such a precise statement. Once Theorem 4.3 is proved, the second part of the analysis [JNRZ14] 
may also be applied to (1.1). This yields a very precise description of the large-time asymptotic 
behavior in terms of a slow modulation in local parameters varying near constant parameters of 
the original wave, and obeying some averaged system of partial differential equations, as derived to 
various order of precision in [Whi74, Ser05, NR13, JNRZ14]. For the sake of conciseness, we do not 
state such a result here but rather refer the reader to [JNRZ14, Theorem 1.12] and accompanying 
discussions in [JNRZ14, Rodl3]. 


5. The shock wave case 

As a sample of the potential wider use of the strategy expounded here, we next turn to the 
connection with viscous shock theory, showing that the same linear damping estimate (3.6) may 
be obtained by essentially the same argument in the asymptotically-constant, viscous shock wave 
case, thus recovering the bounds established in [Zum04, Zum07, GMWZ06] by related but slightly 
different weighted Kawashima-type energy estimates.^ The equations of isentropic gas dynamics in 
Lagrangian coordinates, expressed in a comoving frame are 

dtT — cdxT — dxU = 0 , 

dtu - cdxU dxpir) = ndxiT~^dxu), 

where r is specific volume, u is velocity, and p is pressure. 

Traveling waves {T,u){x,t) = (f,u)(x) satisfy the profile ODE 

—c^T — p(t) -h q = cu'jr, q = constant. 

We note as in the periodic case that c 7 ^ 0, else u,p{t) = constant, yielding r = constant, a trivial 
solution. Assume that the shock is noncharacteristic, i.e., —p'{t±) 7 ^ cP, hence t± are nondegenerate 
equilibria and the shock profile decays exponentially to its endstates as x —>■ ± 00 . 

^The weights used in [Zum04, Zum07, GMWZ06] are effectively 4>i = 4>2 4>3, i4>i)x = —c(^^ -/(^^)^0i, 

C'> 1. 



The linearized equations are 
(5.2) Tt - CTx - Ma; = 0, 


Ut - cUx - (ar)a; = u{t ^Ux) 


where a := p'{t) + Define to be a smooth interpolant between ^^|x=±oo such that 


(5.3) 

for some positive 6. 




,aT^ 


_/(fl^) = 0(e-%l) 


(5.4) 


Taking as before £{U) := J + \ 4 > 2 T^'u^ + <p 3 TUx^, we find again 

^£{U{t)) = J { - + acps) rl - {^(p2) ul^ + (</>! - a(/>2 + :^(/>3) 


+ 0((||u||jy2 + ||'r||//i)(||u||j^i + ||r||i2)). 


Taking 
we thus have 

A 

dt 

(5.5) 


)a;+ = 0, (l)i-a 4>2 + ^ 4>3 = 0, (/>i(0) > 0, and 0 < (?!>2 = constant < 1, 


£{U{t)) = - / 


Tx + {^4>2) u: 
'^Wm + IkllLa)) 


+ O ((||t(||j^2 + ||r||jyq 

< —riW\uxx\\‘L2 + |'rx|li2) + C"(||u||j^2 + ||T||j^i)(||t(||jyi + ||r||j;^2), 

for some positive C' and r]', and thereby the same linear damping estimate as in the periodic- 
coefficient case: 

d 


dt 


£{U{t))<-rj£{U{t)) + C\\U{t)\\i, 


( CXT^ CXl"^ \ 

-/(-)j, hence (/>i and remains 


(5.6) 

for some positive rj and C. 

As in the periodic case, a crucial point is that 

J \ V ' V 

bounded, so that £{U) ~ (||ua ;||^2 + ||'rx||^ 2 ) modulo ||t||^ 2 ) a property following in this case by 
exponential decay, (5.3). 

The above may be recognized as exponentially weighted Kawashima-type estimates similar to 
those used in the study of viscous shock stability in [Zum04, Zum07, GMWZ06], reflecting the 
growing analogy between the periodic and asymptotically-constant cases. Actually, part of the 
recent activity of the authors, jointly with others, focused on dynamics around periodic waves, 
and culminating more or less in [JNRZ14], was motivated by the will to put its analysis on a par 
with classical ones on asymptotically-constant waves. With this respect the present contribution 
that provides analytical tools necessary to consider large-amplitude periodic waves of hyperbolic- 
parabolic systems should be compared with [MZ04], where, motivated by some clever ’’transverse” 
energy estimates of Goodman [Goo91] in the study of small-amplitude stability, the treatment of 
large-amplitude viscous shock waves was first carried out. In the reverse direction the Conjugation 
Lemma of [MZ05] on asymptotically constant-coefficient coordinate transformations from asymp¬ 
totically constant- to constant-coefficient systems may be thought as analogous to the classical 
Floquet Lemma on periodic coordinate transformation of periodic- to constant-coefficient systems 
of equations. 


6. Discussion 

At the linear level, for a general second-order hyperbolic-parabolic principal part Ut + AUx = 
{BUx)x: > 0, a Kawashima-type estimate is on an energy combining £{U) := {Ux, A^Ux) -|- 
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{U,KUx), with the lower-order {U,A^U), where AA is symmetric positive definite and K is skew 
symmetric, chosen, where possible, so that 

( 6 . 1 ) ^{A^B + KA)>Q. 

When A and B are constant, as arising from linearization around a constant state, and the original 
nonlinear system admits a strictly convex entropy, (6.1) may be reduced to a simple-looking con¬ 
dition that is satisfied by most of systems of physical interest; see [Kaw83, LZ97]. As a result, for 
small-amplitude waves, a suitable choice of K may typically be achieved globally with a constant 
K. However, for large-amplitude shocks, this can be done typically only near x —)■ ±oo where A is 
symmetrizable, and one needs to recover coercivity in the near field |x| < C in a different way. 

As exemplified here- for the first time in a periodic context, the key to the treatment of large- 
amplitude waves, is to choose the ”symmetrizer” Aq jointly with the ’’compensator” K so that one 
may use a clever choice for Aq to relax constraints on K and vice versa. In the present case. 



for either of the St. Venant or isentropic compressible Navier-Stokes equations. The issue in the 
latter case is that symmetrizability holds in general only in the limits x —)■ ±oo, in the former that 
it holds only on average, but in any case not pointwise. However, we have seen that energy esti¬ 
mates can be recovered by modulating classical symmetrizers and compensators with appropriate 
asymptotically-constant, or periodic exponential weights. 

Mark that our analysis, while apparently quite robust, leaves widely open the question of de¬ 
termining, for general systems, what kind of notion of symmetrizability on average could lead to 
similar periodic-coefficient high-frequency damping. 


[AeM91] 
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